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Abstract. The Cosmic Background due to the integrated 
radiation from galaxies over the whole life of the Universe 
is reviewed. We find that this background is well con- 
strained by measurements. The total power in the back- 
ground is in the range 60-93 nWm _2 sr~ 1 . The data show 
the existence of a minimum separating the direct stellar 
radiation from the infrared part due to radiation reemitted 
by dust. This reemitted dust radiation is about 1-2.6 time 
the background power in the optical/near-IR thus much 
larger than the same ratio measured locally (30%). The 
far-infrared and submillimeter background is likely to be 
dominated by redshifted infrared galaxies. The long wave- 
length spectrum of the background being significantly flat- 
ter than the spectrum of these galaxies it strongly con- 
strains the far-infrared radiation production rate history 
which must increase by a factor larger than 10 between 
the present time and a redshift 1 and then stays rather 
constant at higher redshift, contrary to the ultraviolet ra- 
diation production rate which decreases rapidly. 

Several models of galaxy evolution have been proposed 
to explain the submillimeter background. In this paper we 
do not propose a new model; we systematically explore the 
allowed range of evolution histories allowed by the data. 
If infrared galaxies are mostly powered by starbursts as 
indicated by recent observations, this infrared production 
history reflects the history of star formation in the Uni- 
verse. 



1. Introduction 

The history of star formation in the Universe is one of 
the key function in physical cosmology. It is closely linked 
to galaxy formation and evolution and controls the sec- 
ond most important contribution to the cosmic electro- 
magnetic background after the Cosmic Microwave Back- 
ground (CMB) generated at the time of recombination 
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at a redshift around 1000. It has been pointed out many 
times in the past 30 years, that measurements of the cos- 
mic background radiated by all galaxies over the history 
of the Universe would be extremely valuable for physical 
cosmology. It would strongly constrain models for galaxy 
formation and evolution (see for example Partridge & Pee- 
bles, 1967). This background is expected to be composed 
of three main components: 

- The stellar radiation in galaxies concentrated in the 
ultra-violet and visible with a redshifted component in 
the near InfraRcd (IR) 

- A fraction of the stellar radiation absorbed by dust ei- 
ther in the galaxies or in the intergalactic medium 

- The radiation from active galactic nuclei (a fraction of 
which is also absorbed by dust and reradiated in the far- 
PR). 

The energy in the first two components is derived from nu- 
cleosynthesis in stars, the last one probably derived from 
gravitational energy of accreted matter onto massive black 
holes. In the last two years the cosmic background at visi- 
ble, IR and submillimeter (submm) wavelengths has been 
finally constrained by very deep source counts and upper 
limits on the diffuse isotropic emission at shorter wave- 
lengths, and measured in the submm range. We review 
the observational situation in Sect. 2. In Sect. 3, we de- 
fine the formalism of the determination of the IR radiation 
production rate history. When the spectrum of the sources 
dominating the background is strongly peaked around a 
wavelength of 80 /im, the radiation production rate as a 
function of redshift can be directly inferred from the spec- 
trum of the cosmic background. In Sect. 4, the sources of 
Cosmic Far-IR BackgRound (CFIRB) and their spectra 
are rewiewed. In the likely hypothesis of this background 
being dominated by IR galaxies (either starburst galax- 
ies or dust enshrouded AGNs), we derive the IR radiation 
production rate as a function of redshift in the Universe 
(Sect. 5). We compare it with other means to measure the 
star formation rate in Sect. 6. 



The study presented here is different from the many 
studies of this type based on a model for Spectral Energy 
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Distribution (SED) and evolution of IR galaxies which, 
after some adjustements, account for the Cosmic Back- 
ground. Considering the diversity of SED which can be 
used and the observational uncertainties it is clear that a 
range of histories of the IR galaxies is possible. The goal of 
this paper and of the method presented here is to explore 
systematically the range of possible solutions by a careful 
analysis of possible SEDs taking properly into account the 
observational uncertainties. The method developed in this 
paper is the only one so far to answer this question. 

2. The Cosmic Background Radiation: observa- 
tions 

We summarise here the observational situation of the cos- 
mic background which is the main observation used in this 
study. The values of the cosmic background at all wave- 
lengths are discussed in Appendix 1 and shown on Fig. 1. 
We can compute the energy contained in the background 
above and below 6 ^m using the best fits presented on 
Fig. 1. We obtain: 

- E(<6 ^m)=E(opt)= 2-4.1 1(T 8 W m~ 2 sr" 1 

- E(>6 /Lim)=E(FIR)= 4.-5.2 1CT 8 W mT 2 sr" 1 

We know that for local galaxies the ratio E (FIR) /E (opt) 
is about 0.3 (Soifer & Neugeubauer, 1991). We measure 
here a higher ratio which is about 1-2.6. This higher ra- 
tio results from two combined effects. First, the redshift 
effect will bring the light from optical to IR and thus will 
give more energy in the IR than for local galaxies. Nev- 
ertheless this effect cannot account for an IR to optical 
ratio much larger on average than the value observed to- 
day. This higher ratio can only be explained by a change 
of properties of galaxies that are making the background 
in the optical and the IR. A very simple comparison sup- 
ports this idea: the energy in the backgound at 15 /J,m and 
in the optical domain is nearly the same. But the back- 
ground at 15 /im in the Hubble Deep Field is made by 
a very small number of objects as compared to the num- 
ber of sources which contribute to the bulk of the energy 
observed in the optical (Aussel et al. 1999). This means 
that the background in the IR is probably dominated by 
few objects with high IR luminosities. Thus these objects 
are not the standard IR counterparts of normal spiral and 
irregular galaxies. 



where ip is the comoving energy production rate per unit 
volume at redshift z (note that here ip is not the spectral 
energy density). 

The backgound spectral intensity (brightness per solid 
angle and frequency interval) I„ is such that: 



(2) 



l+z = 



where we assume a monochromatic spectrum for the 
sources radiating at frequency v =v (1+z) where v is the 
observed frequency. In that case (<5 function spectrum for 
the emitting objects), ip(z) can be simply deduced from 
the background spectrum vl v (see Appendix 2). 

The relations given above can be easily generalised to 
the case where the sources radiate through a broad spec- 
trum: 



4tt J v i =v 



N,L, 



dt 

dz 



dv 



(3) 



with N z the number of sources per Mpc 3 and L v ' the lu- 
minosity of the galaxies (in W/Hz). 

L v > is the average over a large volume of SED of luminous 
IR galaxies which have been shown to dominate the 15 fim 
background at z ~ 0.7 (Aussel ct al 1999). The range of 
SED for IR galaxies are discussed in Sect. 4. The impor- 
tant result of this discussion is that all the far-IR SEDs 
of luminous IR galaxies are not very dependent of the 
energy source (starburst or dust enshrouded AGN) and 
varies slowly with the luminosity (Maffei, 1994; Guider- 
doni 1998). 

The unknown quantity in Eq. (3) becomes N z . The full 
range of allowed functions for the number density of 
sources as a function of z is derived using Monte Carlo 
simulations for a set of IR galaxy spectra which covers the 
range of possible ones and for different cosmological mod- 
els. The number density multiplied by the total luminosity 
of each model galaxy gives the IR luminosity density ip in 
function of z: 



ip(z) =N Z L Vo dv 



(4) 



(p(z) depends on the cosmological model via the ^ . 



3. Formalism of the determination of the radiation 
production rate history in the Universe 

For a Universe following a Robertson- Walker metric with 
total density fio) and a cosmological constant f^A, the dif- 
ferential contribution dE to the background energy den- 
sity (monochromatic or integrated over frequency, per unit 
volume) generated around redshift z during time dt is: 



dE 



(i + z; 



dt 



(1) 



4. Sources of the Cosmic Background and their 
spectra 

To determine <p(z), we need to establish the range of ac- 
ceptable average SEDs (Ly ) of the sources that are mak- 
ing the background. 

1 We can note that p(z) x J| is independent of the cosmolog- 
ical model. We nevertheless choose to use ip(z) rather <p(z) x || 
as p(z) correponds at z=0 to a measured value. 
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Fig. 1. Cosmic Background from the UV to the millimeter wavelength. In the UV domain, the upper limit is from Martin et 
al. (1991) and the value from Armand et al. (1994); the optical and near-IR points are from Pozzetti et al., 1998 (triangles) 
and Bernstein et al., in preparation (squares); The 3.5 and 2.2 pim points are from Dwek & Arendt (1998) and Gorjian et al. 
(2000); Squared upper limits are from Hauser et al. (1998) and crossed upper limits from Biller et al. (1998); the upper limit 
"CAT" is from Barrau (1998) and Barrau et al., in preparation; The 6.5 (Desert, private communication), 12 (Clements et al. 
1999) and 15 /im (Elbaz et al 1999) lower limits come from ISOCAM number counts; the value at 15 /xm (o) is an extrapolation 
of the counts using the Guiderdoni et al. (1998) model. At longer wavelength, we have the 100, 140 and 240 /im Lagache et al. 
2000 (A) and Hauser et al. 1998 (o) DIRBE values, lower limit from Dwek et al. (1998), upper limits from number counts at 60 
(Lonsdale et al. 1990), 175 (Puget et al. 1999) and 850 iim (Barger et al. 1999). Dotted lines are an attempt to draw continuous 
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Table 1. Co- moving luminosity density ip derived from the CFIRB for the three galaxy model luminosities, the four dust 
spectral indexes and the three cosmologies (defined by Qa and flo). Numbers on the rigth side of each row correspond to the 
miminal and maximal value of ip compatible with the CFIRB. The Hubble constant is the same for all cases (h=0.65). Given 
in column 10 is the \ 2 of the CFIRB "fit" derived from the luminosity density compared to the FIRAS and DIRBE measured 
values. 
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Locally Soifer & Neugebauer (1991) have studied the pop- 
ulation of IR galaxies at low redshift using the IRAS all 
sky survey. They have established the integrated IR lu- 
minosity per unit volume in the Universe today, its SED, 
the luminosity function of IR galaxies and the SED as a 
function of the integrated IR luminosity. The IR luminos- 
ity function shows two main components. One is associ- 
ated with normal spiral galaxies which radiate a fraction 
of their energy in the IR (our Galaxy for example radi- 



ates about one third of its luminosity in the IR). These 
galaxies are expected to have a luminosity function sim- 
ilar to the luminosity function of galaxies in the optical. 
However, the IR luminosity function does not show an ex- 
ponential cut off like optical galaxies but display a power 
law behaviour at high luminosities: = L~ 2 . This is also 
shown by more recent compilations of the IR luminosity 
functions (Sanders & Mirabel, 1996). This shows that the 
very luminous IR galaxies cannot be the IR counterpart of 
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Fig. 2. Typical spectra of starburst galaxies (Maffei, 1994) for 
different luminosities (continuous line: 3 10 12 L©, dotted line: 5 
10 11 Lq, dashed line: 7 10 10 L Q ) with a spectral dust emission 
index in the far-IR of 2. 

optical galaxies. In fact they have been shown to be star- 
burst galaxies often associated with merging or interacting 
systems and for which the ratio of IR to optical luminos- 
ity increases with the bolometric luminosity (Sanders & 
Mirabel, 1996; Fang et al. 1998). One can thus conjecture 
that the luminosity function of mergers and interacting 
galaxies is very different from the luminosity function of 
normal galaxies. Using this rough separation, locally the 
relative contributions of starburst to normal galaxies is 
less than 10%. The integrated luminosity of normal galax- 
ies is dominated by L* optical galaxies (L/#* ~ 10 10 L©) 
whereas the integrated luminosity of IR starburst is dom- 
inated by galaxies with Lj# ~10 n L© (see for example 
Sanders & Mirabel, 1996). 

26 of the galaxy detected by ISOCAM in the HDF north 
are identified with galaxies with known redshifts. Aussel 
(1999) has shown that the IR galaxies have luminosities 
greater than 3 10 10 L© at 8.5 fxm (in the rest frame) with 
a median redshift of 0.75, and thus bolometric luminosi- 
ties between 1 and 3 10 11 L©. These galaxies, identified 
mostly as interacting systems or spiral galaxies, are lumi- 
nous IR galaxies undergoing a starburst phase. 
Spectra of IR galaxies have been modeled by Maffei 
(1994), using the observational correlation of the IRAS 
flux ratio 12/60, 25/60 and 60/100 with the IR luminos- 
ity (Soifer & Neugeubauer, 1991). Examples of spectra for 
three different luminosities are shown in Fig. 2. The lu- 
minous IR galaxies are emitting more than 95% of their 
energy in the far-IR. Taking only such kind of galaxies ob- 
viously fails to reproduce the optical Extragalactic Back- 
ground (EB). That is why we concentrate hereafter only 
on the far-IR part of the EB (the so-called Cosmic Far-IR 
Background, CFIRB). 



It is possible however that part of the CFIRB energy 
might be due to dust enshrouded Active Galactic Nuclei 
(AGNs) for which the far-IR SED is very similar to star- 
burst galaxies of similar luminosity. Based on the assump- 
tions that 10% of the mass accreting into balck hole is 
turned into energy and that the black hole masses mea- 
sured in the HDF (Ford et al. 1998) are typical of galaxies, 
the AGN background energy would be in order of 10% of 
that from stars (Eales et al. 1999). These calculations are 
highly uncertain but are supported by the recent work of 
Almaini et al. (1999). The SEDs of high-redshift dust en- 
shrouded AGNs are similar to those observed locally, and 
one can explain 10-20 % of the CFIRB. Because of the 
similar far-IR SEDs, the question of the fraction of the 
CFIRB due to AGNs is not relevant to the determina- 
tion of the IR radiation production history. The estimates 
mentioned above are only important because they indicate 
that the IR radiation production history is likely to reflect 
the star formation history. Locally the average far-IR SED 
adjusted on the IRAS data (Soifer & Neugeubauer, 1991) 
is well represented by the SED of a 8 10 10 L© IR galaxy. 
As we have seen, at z ~ 0.7 the mid IR production is dom- 
inated by galaxies typically 2.5 times more luminous. Fi- 
nally at higher redshift the main indication comes from the 
SCUBA deep surveys and also indicates an average SED 
dominated by ultraluminous IR galaxies. To take into ac- 
count the uncertainty on the average far-IR SED, we use 
for the determination of the IR radiation production his- 
tory, the SEDs from Maffei (1994) for galaxies with lumi- 
nosities 7 10 10 L0, 5 10 11 LQ and 3 10 12 L© which takes 
into account the change in the peak emission wavelength 
with luminosity. For the long wavelength behaviour, we 
allow the dust emissivity index of the model to vary from 
1.3 to 2 which covers generously the uncertainty on this 
parameter for the average SED. 



5. Luminosity density history from Monte Carlo 
simulations 

5.1. The CFIRB spectrum 

We used the CFIRB determination as described in La- 
gache et al. (2000). The CFIRB resulting from an integral 
over a significant redshift range is expected to be a smooth 
function of frequency. Therefore, we use a smooth fit for 
the determination of <p(z). The CFIRB spectrum, can be 
fitted between 200 and 2000 [im by a modified Planck 
function as given by: 

I{v) = 8.80 x \{r b {v/v ) 1A P v {li.QK) (5) 

where ^o=100 cm -1 . The set of parameters (T, t, a) has 
been determined by a x 2 minimization. The spectrum is 
sampled at 10 FIRAS frequency. The uncertainties on each 
sampled frequency are obtained by varying parameters un- 
til x 2 is increased by 10%. In addition to the FIRAS data, 
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Fig. 4. Co-moving luminosity density distribution as derived from the CFIRB for different sets of parameters (with a fixed dust 
spectral index of 1.7). a: for a given cosmological model (f2o=0.3, O.a=0-7) and different luminosities (star: 3 10 12 Lq, diamond: 
5 10 11 Lq, triangle: 7 10 10 Lq). b: given a luminosity of 5 10 11 Lq and three cosmological models with h=0.65 (star: fio=0.3, 
£Ia=0.7, diamond: f2o=0.3, Oa^O, triangle: f2o = l, Qa~0). In each plot, points around a redshift value have been arbitrarly 
shifted to see the uncertainties. 
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Fig. 3. Co-moving luminosity density distribution as derived 
from the CFIRB for all cases listed in Table 1. Also shown is 
the best fit passing through all cases. 

we have used the CFIRB DIRBE determinations at 100, 
140 and 240 fim (Lagache et al. 2000). CFIRB data points 
and uncertainties that are used for the determination of 
(p(z) are shown in Fig. 5. 

5.2. Determination of ip(z) 

We determine the range of functions <p(z) allowed by the 
data for each combination of SED and cosmological model. 
The cases used are given by the combinations of: 



— Three IR galaxy luminosities (3 10 12 L Q , 5 10 11 L Q 
and 7 10 10 L©, see Fig. 2) 

— Four values for the dust spectral index (1.3, 1.5, 1.7 
and 2) 

— Three cosmological models defined by the set of param- 
eters h, n and n A (h=0.65, fi =0.3, fi A =0.7; h=0.65, 
O =0.3, fi A =0, and h=0.65 fi =l, ^A=0) which fix 
dt/dz. 



To establish <p(z) and the associated error bars we use 
Eq. (3). The basic algorithm for finding N z is based on 
Monte Carlo simulations. N z is sampled at a few redshift 
values and linearly interpolated between these values for 
computing the term on the right side of Eq. (3). It has 
been assumed that beyond z=13, N z =0. For each case, the 
solution (minimum % 2 ) has been obtained by exploring a 
wide range of randomly distributed values of N z at each 
of the sampled values of z. Error bars are estimated by 
keeping the computed CFIRB within the uncertainties at 
each sampled frequency (see Fig. 5) and greater than the 
lower CFIRB limit at 850 ^m from Barger et al. (1999). By 
using an iterative method for progressively reducing the 
range of explored values, we reach convergence for each 
case studied with a reasonable number of hits (~ 50000) . 
At z=0, the IR production rate in the Universe is around 
1.65 10 8 h L Q /Mpc 3 (Soifer & Neugeubauer 1991), where 
h is the Hubble constant in units of 100 km/sec/Mpc 3 . 
This gives ip(z = 0)= 10 s L Q /Mpc 3 for h=0.65. 
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Fig. 5. CFIRB models induced from the luminosity density shown in Fig. 4a for the three luminosities (continuous line: 3 10 12 
Lq, dotted line: 5 10 11 L , dashed line: 7 10 10 L ), a cosmology h=0.65, fio = l, ^a=0. Also shown are the observational 
constraints and the CFIRB FIRAS and DIRBE spectrum (stars with error bars) computed as explained in Sect. 5.1. 



5. 3. Results 

The results 2 are summarised in Table 1 for the different 
cosmologies and SEDs. We see that for all the considered 
cases, the \ 2 is very similar. This shows that, as expected, 
there is not a unique solution for the inversion in terms of 
cosmological model and average SEDs. Nevertheless the 
remarkable result is that there is a range of redshifts in 
which all acceptable solutions have the same behaviour as 
can be seen in Fig. 3 where the luminosity density varia- 
tion, together with its uncertainties allowed by the data, 
is shown for each case. On the one hand, as could be ex- 
pected, constraints on ip(z) are very weak below redshift 
1 (no cosmic background values at mid-IR wavelengths) 
and above redshift 4 (very low signal to noise ratio of 
the CFIRB spectrum above 800 /im). On the other hand 
between redshifts 1 and 4, the CFIRB gives strong con- 
straints on the history of the far-IR production rate which 



2 Figures and results can be found on G. Lagache's WEB 
page: http://www.ias.fr/iasnv/people.html 



cannot be established from any of the present source sur- 
veys. The luminosity density is about 10 times higher at 
z=l than at z—0 and it is nearly constant up to redshift 4. 
Because of this rather constant behaviour, the level is only 
weakly dependent on the galaxy model spectrum taken 
with different peak wavelength (see Fig. 4) as a change 
in the peak wavelength mainly shifts the function <p(z) in 
z. Moreover, the luminosity density in the redshift range 
1 to 4 is not affected by changing the far-IR dust spec- 
tral index from 1.3 to 2. On Fig. 5 is shown the CFIRB 
models induced from the luminosity density variation for 
the three luminosities and a fixed cosmology. For the far- 
IR part of the spectrum, the different luminosities give 
automatically good fits . This is not necessarily the case 
below 100 /im where the result is very dependent of the 
galaxy model spectrum. With the present observational 
constraints, the best fit is obtained for a luminosity of 5 
10 11 L©. This is very consistent with typical luminosities 
of IR galaxies that are making the background at 15 /im 
(Aussel, 1999). 
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Fig. 6. SFR derived from UV/Vis/Near-IR observations (dia- 
monds: Lilly et al. 1996; triangles: Conolly et al. 1997; squares: 
Madau et al. 1996 and crosses: Steidel et al. 1999) corrected 
for extinction by Steidel et al. (1999) and SFR derived from 
the CFIRB (continuous line, scaled to Ho=50 km s" 1 Mpc -1 
for consistency) 



Table 2. SFR (in M Q yr" 1 Mpc~ 3 ) deduced from the CFIRB 
and from deep submm SCUBA surveys (for Ho=65 km s _1 
Mpc~ 3 ). 



Publications 


Redshift 


SFR 


This paper 
Barger et al. 1999 
Hughes et al. 1998 


z=1.9 

Kz<3 
2<z<4 


0.21±0.10 
0.25 
>0.14 



6. Cosmological implications 

6.1. The star formation history of the Universe 

Many models describing the evolution of galaxies includ- 
ing their IR and submm emission have been published in 
the recent years. In this section, we discuss only empirical 
determinations of the SFR derived from different observa- 
tions. 

The history of the cosmic Star Formation Rate (SFR) 
can be derived from deep optical surveys assuming that 

(1) the stellar Initial Mass function (IMF) is universal, 

(2) the far-UV light is proportionnal to the SFR and 

(3) extinction is negligible. The presence of dust which 
absorbes most of the UV starlight in starburst galaxies 
makes this last assumption highly questionable. The cor- 
rections needed to account for extinction are rather un- 
certain and there is much controversy about the value of 
this correction. Moreover, the SFR deduced from optical 
surveys can be underestimated if there is a significant pop- 



ulation of objects so obscured that they are not detected 
in these surveys. A direct determination of the fraction of 
the stellar radiation reradiated by dust can be obtained 
from the IR/submm surveys if the dust enshrouded AGNs 
do not dominate. However, so far, the catalogues of faint 
submm sources with reliable redshifts are not large enough 
to reconstruct the history of the SFR (see Lilly et al. 1999 
for a first attempt). We have shown that strong constraints 
can be provided by the CFIRB in the redshift range 1 to 4. 
It is particularly interesting to compare our determination 
with (1) the optically-derived SFR, corrected for extinc- 
tion in the same redshift range and (2) the SFR inferred 
from SCUBA submm surveys. However, it is beyond the 
scope of this paper to discuss the cosmological implica- 
tions in term of star formation. Detail discussions exist 
for example in Dwek et al. (1998), Pei et al. (1999), and 
Madau (1997). 

Fig. 6 shows a compilation of the SFR derived from 
UV/Vis/Near-IR surveys, together with our determi- 
nation of the far-IR radiation production rate history 
(through an interpolation). To compare the star forma- 
tion rate with the luminosity density, we have to use a 
conversion factor. We take M S Q ™-i = 7 riff"L & (Guider- 
doni et al. 1998, with a Salpeter IMF), conversion which 
is in good agreement with that derived from Scoville & 
Young (1983) and Thronson & Telesco (1986). The UV 
luminosity density in Fig. 6 has been corrected for extinc- 
tion by Steidel et al. (1999). We see a very good agreement 
between the UV and far-IR luminosity density confirming 
the need for a very large extinction correction. This is 
a strong indication that the population of galaxies that 
are making the submm EB (300-800 /im) seems to be the 
same as the population detected by Steidel et al. (1999) 
in their surveys of Lyman-break galaxies. This would also 
imply that the population of objects so obscured that they 
are not detected in UV/opt/near-IR surveys cannot con- 
tribute for a large fraction of the luminosity density. 

Submm EB deduced SFR can also be compared to 
SCUBA results. Several groups are now conducting deep 
and ultradeep blanck- field surveys (Barger et al. 1999; 
Hughes et al. 1998; Eales et al. 1999), that follow the 
first survey of Smail et al. (1997) who discover a popu- 
lation of luminous galaxies emitting at 850 /tm amplified 
by lensing from foreground clusters. Two of these groups 
(Hughes et al. 1998; Barger et al. 1999) currently give 
estimates of the submm source SFR whereas Lilly et al. 
(1999) discuss its probable behaviour. Table 2 compares 
the SFR derived from the submm EB and the SCUBA de- 
termination. There is a very good agreement between the 
two although the SCUBA estimates should be interpreted 
with caution since only 20-25 % of the detected sources 
have secure identifications (Sanders, 1999). 
Fig. 5 shows that all galaxies that are contributing to most 
of the background at 850 fim cannot have a luminosity 
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greater than 2-3 10 12 L©, standard luminosities of the 
current detected SCUBA sources (Lilly et al. 1999; Eales 
et al. 1999). This is consistent with the fact that SCUBA 
sources above 3 mJy account for only 20-30 % of the EB 
at 850 /im. The present data show that the bulk of the 
submm EB is likely to reside in sources with 850 /im fluxes 
near 1 mJy. Barger et al. (1999) estimate that the FIR lu- 
minosity of a characteristic 1 mJy source is in the range 
4-5 10 11 L Q , which is what is expected from Fig. 5. 

It can be checked that models of galaxy evolution 
which fit the CFIRB fall within the allowed range of SFR 
histories obtained here (see for example Pei et al., 1999 
and Guiderdoni et al., 1998). 

6.2. Redshift contribution 
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Fig. 7. Relative contribution to the CFIRB of galaxies with 
luminosity of 5 10 11 L for z<2.7 (dotted line), 2.7<z<4.8 
(dashed line), and z>4.8 (continuous line). 

On Fig. 7 is shown the relative contribution from 
different redshift ranges to the CFIRB as a function of 
wavelengths (for one illustrative case a— 2, L=5 10 11 L Q , 
f2n=0.3 and Q\—0.7). As expected, galaxies below redshift 
2.7 contribute mostly at short wavelengths and galaxies 
above redshift 2.7 contribute mostly at long wavelength. 
As can be seen in Fig. 4, the small error bar at z~2.5 
shows that this conclusion is firmly established. On the 
contrary, the fraction shown in Fig. 7 for redshift larger 
than 4.8 is not well constrained by the background. 

At 850 /mi, about 20% of the CFIRB light comes from 
galaxies below redshift 2.7 and 80% from above. This can 
be compared with the recent results of Eales et al. (1999). 
Their data suggest that at least 15% of the 850 /im EB 
is emitted at z<3, which is in good agreement with what 



we obtain here. However, our results show that a signifi- 
cant fraction of the submm EB comes from very distant 
objects. Such a scenario represents a picture in which a 
significant fraction of all stars has been formed very early 
in the Universe. 

7. Summary 

The CFIRB detected in the COBE data at wavelength 
greater than 100 /jm contains a surprisingly large fraction 
of the cosmic background due to distant galaxies. The 
spectrum of this background at long wavelength is sig- 
nificantly flatter than the one observed for individual IR 
galaxies. This implies that the submm part of the CFIRB 
cannot be dominated by the emission of the galaxies which 
account for most of the CFIRB at 150 /im, and thus con- 
tains a unique information about high redshift IR galax- 
ies. Considering the variety of long wavelength spectra 
observed for these galaxies we have explored the range of 
possible redshift evolution histories. We show that only a 
co-moving production rate of far-IR radiation with strong 
evolution at low redshifts but little evolution between red- 
shifts 1 and 4 is the only solution allowed by the CFIRB 
(the detailed low redshift evolution is much better con- 
strained by the ISO deep surveys than by the background) . 
Our results show that there is no evidence for a "peak" 
in the cosmological star formation density at z=l-2 as it 
has been assumed by many authors; it is clear that the 
epoch of the beginning of star formation has not yet been 
identified. Moreover, these results indicate that there is 
a divergence of the behaviour of the star formation his- 
tory as compared to that of the space density of luminous 
AGNs (see for example Dunlop, 1997 and Shaver et al. 
1998), as it was aslo suggested by Steidel et al. (1999). 

Acknowledgements. We thanks Martin Harwit for very useful 
discussions on the content of this paper. The work presented 
here profited from very usefull comments from the referee. 

Appendix 1: Observational constraints on the cos- 
mic background 

At wavelengths at which the contribution is not negligible 
the CMB as well as its dipole distorsion have been removed 
using the values given in Mather et al. (1994) and Fixsen et 
al. (1994). In the following, cosmic background is used to re- 
fer to the extragalactic electromagnetic radiation background 
outside the CMB. 

Far-IR and submillimeter range (A > 100 /im): 

In this wavelength domain, the determination of the cos- 
mic background has been done using the COBE FIRAS and 
DIRBE data. The main difficulty is to remove the bright and 
fluctuating Galactic dust emission. The first detection has been 
reported by Puget et al. (1996). They used an independent 
dataset (the HI 21cm survey of Hartmann et al. 1994) in ad- 
dition to FIRAS spectra to remove the dust emission. The 
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detection of the background at A > 140 /im has been later con- 
firmed by Fixsen et al. (1998), Hauser et al. (1998), Schlegel 
et al. (1998), Lagache et al. (1999) and Lagache et al. (2000), 
using the DIRBE and FIRAS data. All determinations are to- 
day in good agreement for A >150 fim. Discrepancies at 140 
fim are due to different evaluations of the emission from dust 
associated with the ionised gas in the interstellar medium 3 . At 
100 fim, assuming an accurate subtraction of the zodiacal emis- 
sion, Lagache et al. (2000) using two independent gas tracers 
for the HI (Leiden-Dwingeloo survey) and the H + (WHAM 
data) obtain IcF/fls(100)= 0.78±0.21 MJy/sr. One has to 
note that methods based on the intercept of the far-IR/HI 
correlation for the determination of the CFIRB are dangerous. 
For example, for the parts of the sky covered by WHAM, this 
intercept is about 0.91 MJy/sr, which is significantly different 
from the value of the CFIRB (0.78 MJy/sr). It is interesting to 
compare the CFIRB value of 0.78 MJy/sr to the non-isotropic 
residual emission found by Hauser et al. (1998). The average 
over three regions of the residual emission, equal to 0.73±0.20 
MJy/sr, is in very good agreement with the Lagache et al. 
(2000) determination. 

Mid IR range (10 < A < 80 /mi): 

In this wavelength domain, the level of the interplanetary 
dust emission is very high with respect to the expected level 
of the cosmic background. This interplanetary dust emission 
peaks around 25 /im. Around this wavelength, since the level 
of such an emission is around a thousand times the level of the 
background, a direct determination of the background from an 
earth orbit is impossible. At 25 and 60 /tm, our description of 
the interplanetary dust emission is not accurate enough to iso- 
late the background emission. Therefore, in the Mid IR range, 
only indirect determinations of the background are possible. 
Lower limits are provided by deep cosmological surveys. We 
have reported on Fig. 1 lower limits derived from surveys car- 
ried out with the ISO CAM instrument aboard the IR Space 
Observatory: at 6.5 /tm (Desert, private communication), 12 
/mi (Clements et al. 1999) and 15 /im (Elbaz et al. 1999). At 
15 /tm, two values are reported: a lower limit and an absolute 
value. The lower limit is derived from the integrated brightness 
for sources above 50 /iJy. As the counts flattens significantly 
below 100 /tJy it is interesting to display the integrated bright- 
ness down to very low fluxes using an extrapolation of the 
counts with a model fitting the data (Guiderdoni et al. 1998), 
as the result does not depend much on the model used. We 
have also reported a lower limit at 60 /im from IRAS galaxy 
counts (Lonsdale et al. 1990). 

Upper limits are provided by the propagation of very high 
energy gamma rays (in the TeV range) through intergalactic 
space. The principle of such a determination is the following: 
gamma rays are absorbed by pair creation on the photons of 
the cosmic background. For a gamma ray of energy E com- 
ing from a distant extragalactic source all cosmic background 

2 2 

photons with energy above a threshold e = ™g c will con- 

3 The CFIRB value of Lagache et al. (1999) is smaller than 
that derived in Lagache et al. (2000) since the assumed WIM 
(Warm Ionised gas) dust temperature was overestimated (the 
WIM dust spectrum was very noisy below 200 /tm and the 
estimated dust temperature was too high). 



tribute to the absorption of the gamma ray. As more photons 
contribute to the absorption when the gamma ray energy in- 
creases the gamma ray spectrum should present a curvature at 
high energy if the integrated density of photons is high enough. 
Early determinations (Stecker & De Jager, 1993; Biller et al. 
1995) relied on assumptions on the shape of the spectrum at 
the source (a power law extrapolated from lower energies). Re- 
cently variations of the spectrum of the source Mkr 501 have 
been detected with CAT in a burst mode during which the 
X-ray emission has been observed. This allows the prediction, 
with minimal hypothesis, of the initial source spectrum and 
thus get a firm determination of an upper limit (Barrau 1998; 
Barrau et al, in preparation). This limit is in very good agree- 
ment with that of Funk et al. (1998) from HEGRA data. This 
leads to upper limits only a factor of 2 above the lower limit 
given by the extrapolated source counts. 

Ultra-violet, visible and near-IR bands: 

Pozzetti et al. (1998) have compiled number counts of ex- 
tragalactic sources in the Hubble Deep Field (HDF) in the 
U, B, V, I and K bands. The integrated background of these 
counts have converged well and thus we already have a good 
determination of the background due to extragalactic sources. 
Part of the stellar radiation has been absorbed by dust and 
the corresponding energy has been observed in form of the IR 
cosmic background. Another part may have been scattered for 
example by intergalactic hydrogen clouds into a diffuse back- 
ground. Scattering cross sections (by atoms, molecules or dust) 
decrease fast with wavelength thus any scattered isotropic com- 
ponent should be detectable primarily at ultraviolet wave- 
lengths. Direct measurements of the extragalactic background 
are very difficult in the visible range because at these wave- 
lengths, the scattered solar radiation by interplanetary parti- 
cles is larger than the expected cosmic background by a factor 
one hundred. Fortunately in the ultraviolet the fast decrease 
of the solar spectrum allows a good determination of an upper 
limit (Bowyer 1991; Martin et al. 1991). This limit at 1500A 
is shown in Fig. 1 and is at comparable level as the lower limit 
from the HDF counts in U band (at 3600A). Also reported in 
Fig. 1 is the determination of Armand et al. (1994) which is 
in line of the other determinations. The Armand et al. value 
is based on the extrapolation of galaxy counts at 2000A based 
on the balloon-borned experiment FOCA (Milliard et al. 1992). 
We thus conclude that the scattered part cannot exceed 50 % 
in the ultraviolet and is very likely to be negligible in the B, 
V, I and K bands. We can thus make the conservative hy- 
pothesis that the HDF counts provide measurements of the 
cosmic background. However, recent direct measurements of 
the optical background at 3000, 5500 and 8000 A from abso- 
lute surface photometry by Bernstein et al. (in preparation, 
quoted by Madau & Pozzetti, 2000) lie between a factor of 2.5 
to 3 higher than the integrated light from galaxy counts, with 
an uncertainty that is largely due to systematic rather than 
statistical errors. This direct determination is very difficult be- 
cause of the uncertainties level of the zodiacal background to 
be subtracted. 

Also reported in Fig. 1 are the direct estimation of the back- 
ground of Dwek & Arendt (1998) with the DIRBE 3.5 /im 
band and Gorjian et al. (2000) at 2.2 and 3.5 /tm. This de- 
termination is very difficult due to the contamination by the 
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emission of stars. 

The upper limit from high energy gamma rays at 10 /im from 
CAT combined with the extrapolation of the integrated flux 
from the 15 fim counts, and the decrease of the integrated 
emission from sources in the ultra deep ISOCAM surveys be- 
tween 15 and 6 /im, give an indication that the background 
spectrum must present a minimum between 3 and 10 pm. This 
minimum is expected to be the wavelength at which the direct 
stellar contribution becomes equal to the dust (in fact PAHs) 



Appendix 2: Far-IR radiation production rate for a 
single frequency emitting galaxy 

If one makes the extreme assumption that the radiation of 
IR galaxies is concentrated at a single frequency, the spectrum 
of the background can be easily inverted to get the radiation 
production rate. If the production rate of IR radiation per co- 
moving volume element at a single frequency v is <p Vo (z), the 
resulting cosmic background in the simple case of Euclidean 
Universe will be given, following Eq. 2, by: 



L Q Mpc- 



6.510' 



W-»Wm- 2 s 



-5/2 



For a simple approximation of the background with power 
laws 

vl v = 12W- 8 (v/v ) 2 - 5 Wm- 2 sr- 1 
valid between 200 /im and 500 /im, and 

vl v = 2710" 8 (^/i/ o ) 3 Wm" 2 sr _1 

valid between 500 /im and 1.5 mm, (y = 310 12 Hz), we deduce 

4> Vo (z) = 7.8 lO 8 L Mpc~ 3 

for 1 < z < 4 and 

4>„M = 1.710 9 (l + z)-°- B (^^-)°- 5 L Mpc- 3 
for z > 4. 

There is no dependance of the energy production rate with v 
for z between 1 and 4 because the slope of the background 
corresponds to the case of no evolution of the co-moving radia- 
tion production rate. This illustrates why ip(z) does not depend 
very much on the spectrum of the galaxy model in this redshift 
range. 
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